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We study a classically chaotic system which is described by a Hamiltonian Ti.{Q, P; x) where 
(Q, P) are the canonical coordinates of a particle in a 2D well, and a; is a parameter. By changing 
X we can deform the 'shape' of the well. The quantum-eigenstates of the system are \n(x)). We 
analyze numerically how the parametric kernel P{n\m) = |(n(x)|m(xo))P evolves as a function 
of Sx = (x—Xo)- This kernel, regarded as a function of n — m, characterizes the shape of the 
wavefunctions, and it also can be interpreted as the local density of states (LDOS). The kernel 
P{n\m) has a well defined classical limit, and the study addresses the issue of quantum-classical 
correspondence (QCC). Both the perturbative and the non-perturbative regimes are explored. The 
limitations of the random-matrix-theory (RMT) approach are demonstrated. 



I. INTRODUCTION 

Consider a system whose total Hamiltonian is 
T-i{Q, P]x), where {Q,P) is a set of canonical coordi- 
nates, and a; is a constant parameter. This parameter 
may represent the effect of some externally controlled 
field. We assume that both Ho — TLo{QtP\xq) and 
TL = ^.{Q, P; x) generate classically chaotic dynamics of 
similar nature. Moreover, we assume that 6x = (x—xq) is 
classically small, meaning that it is possible to apply lin- 
ear analysis in order to describe how the energy surfaces 
'H{Q, P;x)=E are deformed as a result of changing the 
value of X. Quantum mechanically, we can use a basis 
where Ti.Q = Eo has a diagonal representation, while 

n = Eo + fa B (1) 

For reasonably small h, it follows from general semiclassi- 
cal considerations , that B is a banded matrix. Gener- 
ically, this matrix looks random, as if its off-diagonal el- 
ements were independent random numbers. 

It was the idea of Wigner |^ forty years ago, to study 
a simplified model, where the Hamiltonian is given by 
Eq.(|l]), and where B is a random banded matrix. This 
is known as Wigner's banded random matrix (WBRM) 
model. The applicability of such a model is a matter of 
conjecture. Obviously this conjecture should be tested^. 
The most direct way to test it, which we are going to 
apply, is to take the matrix B of a 'physical' Hamilto- 
nian, and then to randomize the signs of its off-diagonal 
elements. The outcome of such operation will be referred 
to as the effective WBRM model that is associated with 
the physical Hamiltonian. One issue of this paper is to 
make a comparison between the eigenstates of the phys- 
ical Hamiltonian, and those of the associated effective 
WBRM model. 

The standard WBRM model (unlike the 'effective' one) 
involves an additional simplification. Namely, one as- 
sumes that B has a rectangular band profile. The the- 
ory of eigenstates for the standard WBRM model is well 
known Increasing x, starting from Sx = the 

eigenstates of Eq.(|^) change their nature. The general 
questions to address are: 



1. What are the parametric regimes in the parametric 
evolution of the eigenstates; 

2. How the structure of the eigenstates changes as we 
go via the subsequent regimes. 

Recently some ideas have been introduced 1^-01 how to go 
beyond Wigner's theory in case of physical Hamiltonians. 
It has been suggested that there are at least three generic 
parametric scales Sx^"^ <^ Sx^.t <^ 6xsc that control the 
parametric evolution of the eigenstates. We shall define 
these parametric scales later. Accordingly one should 
distinguish between the standard perturbative regime 
{Sx <^ Sx^'"), the core-tail regime {Sx'^"' <^ Sx <$:^ Sx^.t), 
and the semiclassical regime (Sx ^ Sxsc)- 

The purpose of this paper is not just to numerically 
establish (for the first time) the existence of the para- 
metric regimes suggested in l^j-^, but mainly to address 
question (2) above Namely, we would like to study 
how the structure of the eigenstates changes as we go 
via the subsequent regimes. In particular we would like 
to understand the significance of RMT assumptions in 
the general theoretical considerations. The latter issue 
has been left unexplored in the 'quantum chaos' litera- 
ture. (Note however that literally the same question is 
addressed in numerous publication once spectral statis- 
tics of eigenvalues, rather than eigenstate structure, is 
concerned). We also suggest a new procedure for 're- 
gion analysis' of the eigenstate structure. We are going 
to distinguish between first-order tail regions (FOTRs), 
higher-order far-tail regions, and non-perturbative (core) 
region. Our main conclusion is going to be that RMT is 
inadequate for the analysis of any features that go beyond 
first-order perturbation theory. 



II. THE MODEL HAMILTONIAN 

We study the Hamiltonian 

n{Q,P;x) = UP?+Pi + Ql+Ql) + ^ ■ QlQl (2) 
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with X = xq + Sx and Xq = 1. This Hamiltonian de- 
scribes the motion of a particle in a 2D well (see Fig.l). 
The units are chosen such that the mass is equal to one, 
the frequency for small oscillations is one, and for Sx = 
the coefficient of the anharmonic term is also one. The 
energy E is the only dimensionless parameter of the clas- 
sical motion. Our numerical study is focused on an en- 
ergy window around i? ~ 3 where the motion is mainly 
chaotic. 

In the classical analysis there is only one parametric 
scale, which is Sx",} ~ 1. This scale determines the regime 
of (classical) linear analysis. For 6x <C Sx"^ the deforma- 
tion of the energy surface Ho{Q,P;x) = E can be de- 
scribed as a linear process. Later we are going to give a 
precise mathematical formulation of this idea. From now 
on assume that we are in the classical linear regime. 




FIG.l: Left: equipotential contours of the model Hamiltonian 
(^) with X = Xq = 1. Right: A Poincare section of a long 
trajectory (0 < t < 1300) that we have picked in order to 
get the fluctuating quantity J^{t). The initial conditions are 
(Qi,Q2,J'i,^'2) = (1,0,1,2) corresponding to _E = 3. The tra- 
jectory is quite ergodic. It avoids some small quasi-integrable 
islands (the main one is around (0,0)). 

Let US pick a very long ergodic trajectory {Q{t), P{t)) 
that covers densely the energy surface E. See Fig.l. Let 
us define the fluctuating quantity 

^(t) = -{dn/dx) = -Q2g2 

For the later analysis it is important to know the distri- 
bution of the variable J-^ and to characterize its tem- 
poral correlations. The average value is F (JT). 
The angular brackets stand for microcanonical average 
over (Q(0),P(0)), which should be the same as time 
(t) average (due to the assumed ergodicity). The auto- 
correlation function of J-{t) is 

C{t) ^ {{Ht)-F){Ht+T)-F)) (4) 

Note that C(r) is independent of t, and that average over 
t should give the same result as a microcanonical average 
over (Q(0),P(0)). 

The variance of the fluctuations is C(0) = {{J- — F)"^). 
The correlation time will be denoted by r^i. Note that 
with our choice of units ~ 1.0 within the energy range 
of interest. The power spectrum C(u>) of the fluctuating 
J-{t), is obtained via a Fourier transform of C(r). See 
Fig. 2. The average F and the variance C(0) determine 
just the flrst two moments of the J- distribution. The 
probability density of T will be denoted by Pp [J') ■ 



All the required information for the subsequent semi- 
classical analysis is contained in the functions C(r) and 
Pf(^) as defined above. All we have to do in order to 
numerically determine them is to generate one very long 
ergodic trajectory (see Fig.l), to compute the respective 
T{t), and from it to extract the desired information (see 
Fig. 2 and Fig. 3). It is convenient to express Pf{^) in 
terms of a scaling function as follows 

"'<^> ^ vm'^- {-7m) 

By this definition the scaled distribution P^if) is char- 
acterized by a zero average ((/) = 0), a unit variance 
((/^) = 1), and it is properly normalized. Note that 
Pci{—f) rather than Pci(/) correspond to Pf(J^). This 
has been done for later convenience. 



III. THE QUANTIZED HAMILTONIAN 

Upon quantization we have a second dimensionless pa- 
rameter h. For obvious reasons we are considering a de- 
symmetrized (1/8) well with Dirichlet boundary condi- 
tions on the fines Qi=0 and Q2=0 and Qi=Q2- The 
matrix representation of 7i = TL{Q, P;x) in the basis 
which is determined by 7i(Q,P;0) is very simple. The 
eigenstates {n — 1,2,3- ••) of the chaotic Hamiltonian 
TCo = Ti-iQ, P; 1) has been found numerically. 

The phase space volume (dQdP integral) which is en- 
closed by an energy surface 7i{Q, P;x) = E is given 
by a function n = fl{E^x). It is convenient to mea- 
sure phase space volume in units of (27r?i)'', where d = 2 
is the dimensionality of our system. Upon quantization 
the phase space volume n corresponds to the level index 
{n — 1, 2, 3 • • •). This is known as Weyl law. It follows 
that g{E) = dE^{E,x) corresponds to the density of 
states, and A = l/g{E) oc is the mean level spacing. 

In the following presentation we are going to assume 
the our interest is restricted to an energy window which 
is 'classically small' but 'quantum mechanically large'. 
In the numerical analysis of our model Hamiltonian the 
energy window was 2.8 < E < 3.1, where the classical 
motion is predominantly chaotic. The mean level spac- 
ing for ~ 3 is given approximately by the formula 
A w 4.3 * . Our numerical analysis has been carried 
out for h — 0.03 and for h ~ 0.015. Smaller values of h 
where beyond our numerical capabilities since the maxi- 
mal matrix that we can handle is of size 5000 x 5000. 

The representation of in the basis which is de- 

termined by the chaotic Hamiltonian TYq, gives the ma- 
trix B of Eq.([|). The banded matrix B and the band 
profile are illustrated in Fig. 2. The band profile is im- 
plied by the semiclassical relation j^]: 

As we see from Fig. 2 the agreement with this formula 
is remarkable. For the bandwidth Eq.(^) implies that 
Afc ~ 2TTh/T^i. It is common to define 6 = Ah/A. 
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FIG. 2: The band profile {2iTh/ A)-\B„m\^ versus a; = {E„-E,n)/h 
is compared with the classical power spectrum C'{uj). Inset: An 
image of a piece of the B matrix. 



IV. DEFINITION OF THE LDOS PROFILE 

The quantum-eigenstates of the Hamihonian 
TL{Q,P;x) are \n{x)), and the ordered eigen-energies 
are En{x). We are interested in the parametric kernel 

P{n\m) = \{n{x)\m{xo))\'^ = trace(p„p„) (7) 

In the equation above pm{Q,P) and pn{Q,P) are the 
Wigner functions that correspond to the eigenstates 
\m{xo)) and \n{x)) respectively. The trace stands for 
dQdP/ {2TThY integration. 

We can identify P{n\m) as the local density of states 
(LDOS), by regarding it as a function of n, where m is 
considered to be a fixed reference state. An average of 
P{{m + r)\m) over several m-states leads to the LDOS 
profile P{r). Alternatively, fixing n, the vector P{n\m) 
describes the shape of the n-th eigenstate in the Ho repre- 
sentation. By averaging P{n\{n — r)) over few eigenstates 
one obtains the average shape of the eigenstate (ASOE). 
The ASOE is just P(-r). Thus the ASOE and the LDOS 
are given by the same function. One would have to be 
more careful with these definitions if TLq were integrable 
while 7i non-integrable. 

The kernel P(n\m) gives the overlap between the nth 
eigenstate oiH and the mth eigenstate of TYq. For Sx — 
we have simply P{n\'m) — Snm- For dx > the kernel de- 
velops a structure, which is described by the LDOS pro- 
file P{r). If 6x is very small then evidently P{r) consists 
of Kronecker delta (at r = 0) and tail regions (|r| > 0). 
Later we are going to distinguish between first-order tail 
regions (FOTRs), and higher order far-tail regions. As 
Sx becomes larger a non-perturbative core region appears 
around r = 0. Namely, the profile exhibits a bunch of 
states (rather than one) that share most of the probabil- 
ity. If dx becomes even larger, the distinction between 
core and tail regions become meaningless, and the LDOS 
profile becomes purely non-perturbative. We are going 
to explain that the non-perturbative profile reflects the 
underlying classical phase space structure. 



V. THE CLASSICAL APPROXIMATION FOR 
THE LDOS 

The classical approximation for P{n\m) follows 

naturally from the definition Eq.(0). It is obtained if we 
approximate pn{Q,P) by a microcanonical distribution 
that is supported by the energy surface TLiQ, P; x) — 
En{x). Namely, 

Pn{Q,P) - -l—6{'H{Q,P;x)^E,,(x)) 

= 5inin{Q,P;x))-n) (8) 

and a similar expression (with x = xq) for PmiQ, P)- In 
the classical limit n is the phase space volume by which 
we label energy surfaces. Each energy surface n is associ- 
ated with a microcanonical state Pn{Q, P)- The classical 
LDOS profile will be denoted by P^i{r). The Sx regime 
where the classical approximation P{r) « Pci('') applies 
will be discussed in a later section. 

By definition, for Sx <^ Sx"^ the deformed energy sur- 
faces departs linearly from the dx = surfaces. As al- 
ready stated in the Introduction, being in this classical 
linear regime is a fixed assumption of this paper. Now 
we want to explain the consequences of this assumption. 
One may consider these consequences as giving an op- 
erational definition for the classical linear regime. The 
dispersion (square-root of the variance) of the classical 
profile in the classical linear regime is 

SE,, = VC{0j X 6x (9) 

(This should be divided by A if we want the disper- 
sion in proper r units. See ( pd] ) below). For our model 
Hamiltonian, for energies £^ ~ 3, we have found that 
SE^i « 0.38 * Sx. Eq.(|^) can be regarded as a special con- 
sequence of the following scaling relation which we are 
going to derive below: 

PAr) - -^^■PJ^^=^] (10) 
,/C{0) Sx \ VC(0) SxJ 

The scaling function has already been defined in Eq.(ra), 
and it is illustrated in Fig. 3. The classical profile Pdir) 
is in general non-symmetric, but it follows from Eq.(p^) 
that it must be characterized by (r) = 0. [By defini- 
tion the scaling function of Eq.(||) gives zero average]. 
Another obvious feature is having sharp cutoffs, beyond 
which Pci{r) = 0. The existence of these outer 'classically 
forbidden' regions follows from the observation that for 
large enough r there is no longer classical overlap between 
the energy surfaces that correspond to \m{xo)) and \n{x)) 
respectively. 

The rest of this section is dedicated to technical clarifi- 
cations of Eq. ([lO|) , and it can be skipped in first reading. 
The derivation is done in two steps. The first step is to 
establish a relation between P^iir) and its trivially re- 
lated version Psii)- The second step is to demonstrate 
that PeIs) is related to Pf{^) of Eq.(|^). It is also pos- 
sible to make a one-step derivation that relates Pci('') to 
Py{J-), but we find the derivation below more physically 
appealing. 
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FIG. 3: The scaled classical profile PdO- One unit on the hori- 
zontal axis correspond to energy difference SE^i ^ 0.38 * 5x. Note 
that r = implies (En{x)—Em{xo)) > 0. The caustic is lo- 
cated at (_E„(x)— _Em(a;o)) = 0, while the anti-caustic is located at 
{En(x)—Em(xo)) = 1.65 * X. The "forbidden regions" are defined 
as those regions where Pci{r) = 0. They are located to the left of 
the caustic and to the right of the anti-caustic. 

By differentiating of n = Vl{E,x)^ keeping n con- 
stant, we get the relation SE = —F{x)5x^ where F(x) = 
dx^{E,x)/g{E) is known as the (generahzed) conserva- 
tive force. Using the latter expression it is a straightfor- 
ward exercise to prove that F{x) — {J-) = F. Alterna- 
tively, we can eliminate E from the relation n = x), 
and write the result a.s E = En{x). Accordingly F{x) = 
— {dEn{x)/dx). Now we can write the following relation: 



En{x) - Emixo) 



dE 
dx 



Sx 



dE 
dn 



(n — to) 



which can be re-written in the following form 

e = -F{x) Sx + {l/g{E)) r 



(11) 



Whenever we regard the kernel P{n\'m) as a function of 
n — m 'we use the notation P{r). But sometimes it is 
convenient to regard P(7i|to) as an energy distribution 
PeCc)- Due to the change of variables ( [ll|) we have the 
following relation: 



P(r) 



1 



P. 



9{E) ^\9{E) 



1 



-r — F{x)5x 



(12) 



The energy distribution P-E{e) can be formally defined as 
follows: 

P^{e) - ^P(n|TO) <5(e-(£;„(x)-£;„(a;o))) (13) 



VI. NUMERICAL DETERMINATION OF LDOS 
PROFILES 

Given 5x we can determine numerically the LDOS pro- 
file P{r). Representative profiles are displayed in Fig. 4. 
For the purpose of further discussion we introduce the 
following definitions: 

• The classical LDOS profile F,,(r) 

• The quantum mechanical LDOS profile P{r) 

• The effective WBRM LDOS profile P^m-rir) 

• The first-order perturbative profile Pp^ {r) 

We have already discussed the classical LDOS profile. 
Below we explain how we numerically determine the 
quantum mechanical LDOS profiles P{r) and Prmt(?'), 
and we also define the profile Pprt(r). 

The numerical procedure for finding P(r) is straight- 
forward. For a given 6x we have to diagonalize the ma- 
trix (^. The columns of the diagonalization matrix Tmn 
are the eigenstates of the Hamiltonian, and by defini- 
tion we have P(ri|m) = |Tm„p. Then P{r) is com- 
puted by averaging over roughly 300 reference states that 
are located within the classically-small energy window 
2.8 < E < 3.1. Fig. 4 displays typical profiles. 

The effective WBRM Hamiltonian is obtained by ran- 
domizing the signs of the off-diagonal elements in the B 
matrix. For the effective WBRM Hamiltonian exactly 
the same procedure (as for P(r)) is applied leading to 

Prmt(^)- 

In order to analyze the structure of either P{r) or 
Prmt (r) we have defined the first-order perturbative pro- 
file as follows: 



r2 + [En-Er, 



(14) 



It is implicit in this definition that (En— Em) and |B„mp 
should be regarded as a function of r. The r = value 
of the band-profile should be re-defined by an interpola- 
tion. The parameter F = 6oA is determined (for a given 
6x) such that the Pprt(f) has a unit normalization. Note 
that Wigner's Lorentzian would be obtained if the band 
profile were flat. 



In the classical limit the summation over n should be 
interpreted as a dn integral. For P(n|TO) in the above 
expression we can substitute the definition Eq.(^ with 
Pn and pm approximated as in Eq.(^). A straightforward 
manipulation leads to the result: 

P,(e) = (Sie- iH{Q,P;x) -HiQ,P;xom 



= {6{e + SxT{t))) = — Pp 




Together with (||) and (|T^ we get Eq.(|T^) along with 
the implied special result(^. 



VII. REGION ANALYSIS FOR THE QUANTAL 
LDOS 

By comparing P(r) to Pp^i{r) as in Fig. 4, we can 
determine''' the range bi [left] < r < bi [right] where Pprt (r) 
is a reasonable approximation for P{r). Loosely speaking 
(avoiding the distinction between the 'left' and the 'right' 
sides of the profile) we shall say that Pp^t (r) is a reason- 
able approximation for |r| < bi. The core is defined as 
the region |r| < 6o- The FOTRs are 6o < |r| < 6i. The 
far-tail regions are |r| > bi. 
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FIG. 4: The quantal profile P{r) is compared with Pp,t(r) and 
with Prmt(''). We are using here the h = 0.015 output. The 
insets are normal plots while the main figures are semilog plots. 
In the lower plot {Sx = 0.2123) the classical LDOS profile Pc]{r) is 
represented by heavy dashed line. 



The results of this region analysis are summarized by 
Fig. 5. In the following sections we are going to present a 
detailed discussion of this analysis. For the convenience 
of the reader we summarize: 

• bo = border of the core region 

• bi ~ border of the first order tail region (FOTR) 

Having 6o 1 implies a standard perturbative struc- 
ture. Having 1 <C &o ^ implies that we have a well 
developed core-tail structure. Having bo ~ bi implies a 
purely non-perturbative structure. In the latter case the 
distinction between core and tail regions become mean- 
ingless. 




VIII. THE STANDARD PERTURBATIVE 
REGIME 

The standard perturbative regime Sx Sx"^"" is de- 
fined by the requirement boldx) 1. This condition 
implies that P{n\m) ~ Snm- For numerical purpose it 
is convenient to define 6x'^"' as the value of 6x for which 
P{r = 0) « 0.5. The theoretical considerations of ||] 
imply that fe^" oc The prefactor is a classical 

quantity whose precise value depends on the operational 
definition of dx"^™. With the operational definition given 
above we have extracted the result fe^" w 3.8 * h"^^^. 

In the standard perturbative regime we can write 
schematically 



P(n\m) 



Tail 



(15) 



The 'Tail' is composed of FOTRs and far-tail regions. 
The former are given by Eq. (|l^) , while the latter are de- 
termined by higher orders of perturbation theory. Note 
that for the standard WBRM we have by construction 
bi = b, and more generally n-th order perturbation the- 
ory becomes essential for (n—l) x b < \r\ < n x b. In 
case of our physical Hamiltonian, as well as for the as- 
sociated effective WBRM model, the boundary bi is 6x 
dependent. 

By comparing P{r) with Prmt('') we can see that RMT 
cannot be trusted for the analysis of the far-tails, because 
system-specific interference phenomena becomes impor- 
tant there. Namely, the RMT profile i-'nMT('') is almost 
indistinguishable from Pprt(r). In contrast to that, the 
far-tails of P{r) are dominated by either destructive in- 
terference (left tail) , or by constructive interference (right 
tail). 



IX. THE CORE- TAIL REGIME 

The core-tail regime Sx'^"' <^ Sx <$:^ Sx^^^ is defined by 
the requirement 1 <C 6o ^ ^i- The theoretical consid- 
erations of 1^ imply that Sx^^t oc h. The prefactor is a 
classical quantity whose precise value depends on the op- 
erational definition of Sx^^-t . In our numerical analysis we 
have defined SXp^t as the 5x for which the contribution of 
the FOTRs to the variance becomes less than 80%. With 
this operational definition we have extracted (using the 
lower subplot of Fig. 5) the result Sxp^t « 5.3 * ?i. 

In the core-tail regime we can write schematically 

P{n\m) w Core -I- Tail (16) 

Disregarding the far-tail regions, the large-scale behavior 
of P{r) can be approximated by that of Pp^t{r). As in 
the standard perturbative regime one observes that the 
far-tails are dominated by either destructive interference 
(left tail), or by constructive interference (right tail). 

The core is a non-perturbative region. It means, that 
unlike the far-tail, it cannot be obtained from any finite- 
order perturbation theory. Once the core appears, the 
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validity of first-order perturbation theory becomes a non- 
trivial matter. In a non-rigorous argument is sug- 
gested in order to support the claim that, disregarding 
smoothing effect, the local mixing of neighboring levels 
does not affect the growth of the tail. An important 
ingredient in this argumentation is the (self consistent) 
assumption that most of the probability is well-contained 
in the core region. Indeed the analysis which is presented 
in Fig. 5 is in agreement with this assumption. 

The observation that the local mixing of neighboring 
levels does not affect the growth of the tail, implies that 
the tail grows as 5x^ and not like say 5x. (The latter 
type of dependence is implied by an over-simplified ar- 
gumentation). Having indeed Sx^ behavior is implied by 
observing that Pir) « P^,t{r) for the FOTRs. 

Finally, it should be emphasized that the local mixing 
of levels on the small scale is not reflected by Eq. ( 14 ) . 
In particular, one should not expect Eq.(14) to be liter- 
ally valid within the core region (jrl < bo). 



X. THE NON-PERTURBATIVE REGIME 



In our numerical study we could not make h small 
enough such that Sx^c <^ ^^t- Therefore, the lower 
profile in Fig. 4 is neither reasonably approximated by 
Pprt(?'), nor by P^i(r). However, we have verified (by 
comparing the h = 0.03 output to the h — 0.015 output) 
that detailed QCC between P{r) and Pd{r) is easily im- 
proved by making Ti smaller. Comparing P{r) to Pci{t) 
on the one hand, and P„i,t(^) to ^'ci('') on the other hand, 
leaves no doubt regarding the manifestation of underly- 
ing classical structures. 

Using a phase-space picture [^,^ it is evident that 
larger 5x leads to better QCC. The WBRM model does 
not have a classical limit, and one finds a quite different 
scenario |^ . For large enough 5x the eigenstates of Eq. (|^) 
become Anderson localized. This localization shows up 
in the ASOE provided the eigenstates are properly cen- 
tered prior to averaging. In the (non-averaged) LDOS 
localization manifests itself as sparsity, and therefore the 
various moments of the LDOS profile are not affected. 
This latter remark should be kept in mind while reading 
the next section. 



In the non-perturbative regime {5x ^ feprt) one may 
say that the core spills over the FOTRs and therefore 
P[n\m) becomes purely non-perturbative. As an exam- 
ple for non-perturbative profile let us consider the lower 
plot of Fig. 4, corresponding to 5x = 0.2123. We see 
that there is poor resemblance between P{r) and Pp,., (r). 
The LDOS profile P{r) no longer contains a predominant 
FOTRs. This claim can be quantified using the analysis 
in Fig. 5. The lower subfigure there displays the FOTR 
contribution to the dispersion. For 5x > Sx^^^ the disper- 
sion is no longer determined by the FOTR contribution. 

The complete disappearance of FOTRs is guaranteed 
only for 5x ^ feprt- Evidently, for 5x :s> Sxp^t the 
FOTRs must disappear, because P{r) goes on expand- 
ing, while Pprt{r) saturates. This is not captured by our 
numerics since for h ~ 0.015 we cannot satisfy the strong 
inequality Sx 3> Sxj,^t, and have a classically small Sx at 
the same time. 



XI. THE SEMICLASSICAL REGIME 

Looking back at the lower plot of Fig. 4, we see that 
detailed QCC with the classical profile (represented by 
heavy dashed line) starts to develop. The right far tail 
contains a component where P{r) and Poi{r) are indistin- 
guishable. This detailed QCC obviously does not hold for 
the RMT profile. 

Being in the non-perturbative regime does not imply 
detailed QCC Detailed QCC means that P{r) 

can be approximated by Pci{r). Having Sx ^ fep,t is 
a necessary rather than sufficient condition for detailed 
QCC. 

A sufficient condition for detailed QCC is Sx ^ Sxsc- 
The parametric scale Sxsc is defined in j^, and for our 
system we can obtain the (theoretical) rough estimate 

Sxsc « 4* 
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FIG. 5 The results of region analysis. The common horizontal axis 
is Sx. The upper subfigure presents the r boundaries as a function 
of 5x. The dotted lines ±60 define the core region (|r| < bo). The 
solid lines define the r region in which 50% of the probability is con- 
centrate. The dashed lines are bi[left] and 61 [right]. The FOTRs 
are the regions where bo < \r\ < bi. The light solid lines and the 
light dashed lines are for the effective WBRM model. The lower 
subfigure displays the dependence of SEci and SE^m and SEpn on 
Sx. The quantal and the classical results are almost indistinguish- 
able, whereas SEprt approaches saturation. The contribution of 
the FOTRs to i5-Eq„ is also displayed. 



XII. RESTRICTED QCC 

It is important to distinguish between detailed QCC 
and restricted QCC. Let us denote the dispersion of the 
quantal LDOS profile by 5E^^, The corresponding clas- 
sical quantity is given by Eq.(^). The two types of QCC 
are defined as follows: 
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Detailed QCC means P{ 



• Restricted QCC means 6E^,„ « SE^i 

Obviously restricted QCC is a trivial consequence of de- 
tailed QCC, but the converse is not true. It turns out that 
restricted QCC is much more robust than detailed QCC. 
In Fig. 5 we see that the dispersion SE^^ of either P(r) 
or Prmt('') is almost indistinguishable from SE^i. This is 
quite remarkable becuase the corresponding LDOS pro- 
files (quantal versus classical) are very different! 

It is important to realize that restricted QCC is implied 
by first order perturbation theory. If we use Eq. and 
take into accound the FOTR dominace which is implied 
by Sx <IC Jxprt, then we get simply 



SE^^ — 



^P(n|m) iEn-E,r^y = 6x^ 



E 



(17) 



where prime indicates omission of the n — m term. Using 
Eq.(^ one realizes that this result is in complete agree- 
ment with Eq.(^. In contrast to that higher moments of 
the perturbative profile are vanishingly small compared 
with the corresponding classical result. The latter fact is 
just a reflection of the absence of detailed QCC. 

One may wonder what happens with Eq. ( [17| ) if we try 
to do a better work, taking into account the core width, 
as well as higher order far-tails contributions. One may 
think that Eq.([l^) is only the lowest order approxima- 
tion, which would imply that restricted QCC should be- 
come worse as Sx grows. However, the latter speculation 
turns out to be wrong. 

We already saw that restricted QCC is implied on the 
one hand (for small Sx) by first-order perturbation the- 
ory, and on the other hand (for large 5x) by detailed 
QCC. Now we would like to argue that restricted QCC 
holds in general. It simply follows from the observation 
that 6E^^ is determined just by the band profile. The 
prove is very simple The variance of P{n\m) is de- 
termined by the first two moments of the Hamiltonian in 
the unperturbed basis. Namely 

SEl^^ = {m\H^\m) - {m\n\mf 
= 6x'^{{m\B'^\m) - (TO|B|m}^) 

Thus, we get the same result as in first order pertur- 
bation theory without invoking any special assumptions 
regarding the nature of the profile. Having SE^„^ that is 
determined only by the band profile, is the reason for de- 
tailed QCC, and is also the reason why restricted QCC 
is not sensitive to the RMT assumption. 
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